Many-Body Entanglement: a New Application of the Full Counting Statistics 
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Entanglement entropy is a measure of quantum correlations between separate parts of a many- 
body system, which plays an important role in many areas of physics. Here we review recent work 
in which a relation between this quantity and the Full Counting Statistics description of electron 
transport was established for noninteracting fermion systems. Using this relation, which is of a 
completely general character, we discuss how the entanglement entropy can be directly measured 
by detecting current fluctuations in a driven quantum system such as quantum point contact. 



I. INTRODUCTION 

Density matrices as a tool for describing quantum- 
mechanical systems when only partial information about 
the quantum state is available were introduced in quan- 
tum mechanics by Landau in 1927 [1]. In this work, 
called "The problem of damping in wave mechanics," he 
was interested in irreversibility of certain quantum me- 
chanical processes, such as the spontaneous decay of ex- 
cited atomic states. Such irreversibility is not inherent 
to quantum mechanics, it arises from a fully reversible 
quantum evolution of a larger system, including the vari- 
ables describing radiation. More generally, Landau was 
concerned with the situation when some variables needed 
to completely describe the system cannot be measured. 
In such cases the imprecision of our knowledge renders 
the quantum state vector ip a useless quantity. Instead, 
a density matrix must be used, which is defined in the 
subspace of the system's Hilbert space spanned by those 
states which can be measured, p = • Pij\i)(j\- 

Another generalization of quantum mechanics in which 
density matrices feature prominently is quantum statis- 
tical mechanics. It was developed, also in 1927, by von 
Neumann 0, as a way to introduce statistical descrip- 
tion in quantum theory. In this approach, a quantum 
system can occupy states forming an orthonormal 
set, with statistical probabilities < pi < 1. Such a sys- 
tem is described by p = $^vPi|£)(i|, which is nothing but 
a diagonal representation of the density matrix. In this 
work, von Neumann first wrote his famous formula for 
the entropy, 



S = -Tr(plo g/9 ) 



Pi log Pi, 



(1) 



which is a proper extension of the Gibbs entropy (and the 
Shannon entropy) to the quantum case. For a set of N 
states, the entropy (JT|) can take values varying between 
(for a pure state) and log TV, realized when all states are 
occupied with equal probabilities. 

The two aspects of the density matrix, emphasized by 
Landau and von Neumann, the lack of information about 
the quantum state in the situation when some variables 
are not measureable, and the connection with statistical 
description, are combined in the notion of entanglement 
entropy This quantity was introduced in 1986 when 



Bombelli et al. [3j] proposed to use von Neumann en- 
tropy, mutual for different parts of a quantum system, 
as a model of the Bckcnstein-Hawking entropy of black 
holes. They considered a situation when a quantum- 
mechanical system is described by a pure state which 
is delocalized between two regions A and B of its con- 
figuration space, such that only one of those regions, A, 
is accessible for physical measurement (see FigO}. An 
observer situated in region A, after performing a full set 
of measurments available to him or her, will describe the 
system by a reduced density matrix pa = T^bP, where 
p = |0)(0| describes the pure state of the whole system 
and Ttb denotes a trace taken over all variables in re- 
gion B. The von Neumann entropy associated with the 
projected state, defined as 



Sa = -Tr PA log pa, 



(2) 



is a characteristic of entanglement between quantum vari- 
ables in the regions A and B. One can easily check that 
entanglement entropy does not change when regions A 




FIG. 1: Entanglement between two regions, A and B, probed 
by an observer who can perform measurements only in region 
A but not in region B. The quantum state of the entire 
system A + B, described by the state vector |0), must be 
"projected" on A, giving pA = TtbPo, where po = |0)(0|. 
Due to entanglement between variables in A and in B, the 
projected state is typically a mixed state, pA ^ Pa, even 
when the whole system is in a pure state. The von Neumann 
entropy of the projected state, Eq.((2}, is a measure of the 
entanglement between A and B. 
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and B are interchanged, Sb = Sa- 

These ideas were developed further by Callan and 
Wilczek [ij], and Holzhey, Larsen and Wilczek [B[ who 
considered a 1 + 1 dimensional system described by con- 
formal field theory. Taking region A to be an interval of 
length £, they found that the entanglment entropy obeys 
the relation 



S = 



G 



log- 



(3) 



where c (c) is a central charge of the conformal field the- 
ory, and a is a microscopic cutoff length. This finding, 
which showed that entanglement entropy is sensitive to 
fundamental characteristics of a quantum system, moti- 
vated many further studies in which entanglement en- 
tropy has been used as a tool to analyze many- body 
states of a variety of different systems [6.7,8,9, lOl UlllTa . 

fUHliHlilliTt. 

Entanglement entropy, serving as a general charac- 
teristic describing quantum many-body correlations be- 
tween two parts of a quantum system, provided a frame- 
work for analyzing quantum critical phenomena 
and quantum quenches [9l fl0lfll|l^ | . Recently it was 
used as a probe of complexity of topologically ordered 
states [ilil6lll~7l | . In addition, this quantity is of fun- 
damental interest for quantum information theory as a 
measure of the resources available for quantum compu- 
tation [l3| as well as for numerical approaches to strongly 
correlated systems [14j . 



II. MEASURING THE MANY-BODY 
ENTANGLEMENT 

The universal appeal of entanglement entropy arises, 
at least in part, from the fact that this quantity is de- 
fined solely in terms of the many-body density matrix 
of the system, with no relation to any particular observ- 
ables whatsoever. This is the main reason this quantity 
has found applications in such diverse fields as cosmol- 
ogy, field theory, condensed matter theory, and quantum 
information. However, for the very same reason, it has 
not been clear how to access this quantity experimentally, 
since measuring the entire density matrix of a many-body 
systems represents a formidable task. Indeed, the many- 
body density matrix depends on coordinates of all parti- 
cles in the system, which usually cannot be measured all 
at once. 

Based on what was just said, the very idea of mea- 
suring a quantity that encodes information about many- 
body correlations of a large number of particles, which 
is what the entanglement entropy is, may seem totally 
bizarre. Yet, the situation with the entanglement entropy 
is different from that of the many-body density matrix. 
Recently we have shown that a direct measurement of en- 
tanglement entropy is possibe owing to its relation with 
quantum noise of electric current KL SI) - 

The system analyzed in Ref.[18j| is a quantum point 
contact, representing an electron beam-splitter with 



transmision and reflection coefficients tunable by exter- 
nal gates [19]. In essence, the QPC serves as a door 
between electron reservoirs, which can be opened and 
closed on demand. By adjusting the voltage on the gates 
of the QPC, different transmission channels in it can be 
opened or closed individually, leading to quantized steps 
in the QPC conductance (see FigJ^-). 

The simplest protocol of driving the QPC that leads 
to entanglement of many-body states in the reservoirs 
is illustrated in FigJSJs. We start with the QPC in a 
closed state, when the resrvoirs are disconnected, then 
the QPC is opened during time to < t < t\, and then 
closed again. In this process, most of the particles either 
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FIG. 2: (a) Quantum point contact (QPC), an electron 
beam-splitter with tunable transmission and reflection (from 
Ref.[l9j]). Conductance G of a QPC exhibits quantized steps, 
observed by varying the gate voltage V g . As V g becomes more 
negative, the transport channels in the QPC open up one 
by one, with transmission Di in channel i increasing from 
to 1 between consecutive steps i and i + 1. (b) Schematic 
of two Fermi seas which are connected via a QPC during 
the time interval to < t < ti, and then disconnected (from 
Ref.[l8j]). Electron transport, taking place at to < t < ti 
makes electrons delocalized among the two leads, generating 
entanglement and current fluctuations. There is a "space- 
time" duality between this situation and the conventional ap- 
proach [KUil, in which many-body correlations are analyzed 
using a finite region in space. 
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remain in their initial resrvoir, or are fully transmitted to 
another reservoir. Some particles, however, become de- 
localized between the two reservoirs, making the states 
in the two reservoirs entangled. After the QPC is closed 
any communication between the reservoirs becomes im- 
possible. This situation mimics that considered in the 
definition of many-body entaglement, when an observer 
can perform measurements only in one region but not in 
the other (see Fig[T]). 

The state of two Fermi seas, coupled via the QPC as 
shown in FigJ^b, evolves as a pure state until projection 
on a specific reservoir is performed. After projecting on 
reservoir L at t = t±, the density matrix takes the form 



pUh) = Tr^U^ioVoUt^o,^)), 



(4) 



where po is the initial state, U is the many-body evolu- 
tion between to and t%, and Tt\r is a partial trace over 
degrees of freedom in the lead R. Due to the exchange 
of particles between reservoirs during to < t < t\, the 
resulting density matrix describes a mixed state with a 
nonzero von Neumann entropy 



St 



-Trp L (*i)logp L (*i) 



(5) 



which characterizes entanglement buildup due to particle 
exchange between reservoirs. 

The evolution of two coupled Fermi seas, describing 
this process, was analyzed in Ref. 18]. This can be done 
exactly owing to the free fermion nature of the problem. 
It was found that all multi-particle correlations in the 
Fermi sea that are relevant for entanglement are fully 
accounted for by temporal correlations of electric current 
flowing through the QPC. Specifically, there is a universal 
relation between the entanglement entropy ([5]) and the 
full counting statistics of the charge transmitted through 
the QPC. This relation, which we shall discuss below, is 
of a completely general nature, independent of the details 
of the protocol used to drive the QPC. As such it can be 
used to obtain the entanglement entropy from measured 
fluctuations of electric current. 

The relation between entanglement and electric noise 
has been at the center of the discussion of different ways 
to generate entangled pairs in a driven electron system, 



using tr ansp ort in normal metal-superconductor junc- 
UHf and in the QPC [HH. Such pairs, which 



tions 

represent an electron analog of the recently demonstarted 
entangled photon pairs [24[ , could be used for testing Bell 
inequalities in a condensed matter system. 

In contrast to Refs. 20.21.22.23], here we are concerned 
with entanglement of many-body states, represented by 
Fermi seas in the right and left reservoirs shown in FigJ^Ja. 
This entanglement is generated by the evolution of the 
full many-body fermion state happening when the QPC 
is opened and then closed. In that, there is an analogy 
with recent literature in which generation of entangle- 
ment in time for critical Hamiltonians d and for generic 
Hamiltonians I(ilfl2| was discussed. 

The centerpiece of the approach of Ref.[l8| is the re- 
lation between many-body entanglement and a physical 



measurement, which in this case is electric current fluc- 
tuations. A relation of entanglement with another mea- 
surable quantity, particle number statistics, was empha- 
sized in Ref. [331] . In this paper a Fermi system with a 
fixed total number of particles was considered in a setup 
pictured in Fig[T] The many-body entanglement in this 
system was expressed through the probability distribu- 
tion of particle number in the region A. Similarly, in the 
approach discussed here, we link the entanglement gen- 
erated in a driven fermion system, such the QPC, to the 
statistics of charge transmitted between reservoirs. 

The relation between entanglement and counting 
statistics of charge fluctuations, which we discuss below, 
was derived in Ref. [HI for a noninteracting fermion sys- 
tem. In the derivation we focus on the QPC as a con- 
venient model, however it will be clear that the result 
is more general. The overall simplicity of the relation 
between entanglement and counting statistics, and also 
its independence of the details of the driving protocol, 
suggests an even higher degree of generality. It would 
be extremely interesting to find out whether a similar re- 
lation holds for interacting many-body systems, such as 
quantum spin chains, Luttinger liquids or Quantum Hall 
liquids. 



III. A PRIMER ON COUNTING STATISTICS 

The Full Counting Statistics (FCS) approach has been 
developed in the theory of quantum noise to describe 
current fluctuations in nanodevices such as QPC and 
tunnel junctions [25[. These fluctuations can be char- 
acterized by the probability distribution of charge trans- 
mitted through the device during the measurment. It is 
convenient to combine individual probabilities in a single 
quantity, the generating function 



X(A) 



E 

n— — oo 



P n e 



iXn 



(6) 



where P n is the probability to transmit n charges. The 
auxiliary variable A is sometimes called "counting field" 
in the literature. 

For example, a binomial distribution with the number 
of attempts N and the probabilities to succeed and fail 
in each attempt p and q = 1 — p is described by 



Pn = 



X(\) = (l-p + pe ix ) 



(7) 



Probability distribution of this form arises in the problem 
of a DC-biased QPC [H. 

The function xW encodes all cumulants of FCS (or, 
irreducible moments) via an expansion 



logx(A) = J2 



(iX) m Cr, 



m l 



(8) 
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The lowest cumulants C\, C2, C3... describe properties 
of the distribution P n such as the mean n, the variance 
((n — n) 2 ), the skewness ((n — n) 3 ), etc. 

The 2nd cumulant Ci is available from routine noise 
measurement. Recently, the 3rd cumulant C3 has been 
measured in tunnel junctions (26l[27l| and in QPC [HI, 
while cumulants up to 5th order where measured in quan- 
tum dots [HHj. In fact, the method used in Refs. [27||2j| 
yields the full probability distribution P n (see Fig(3j) ; 
however, only the lowest moments C\, Ci and C3 of this 
distribution were found to be dominated by intrinsic ef- 
fects. 



N 
H 



< 



C/3 




FIG. 3: Measurement of high-order cumulants of electron 
shot noise in a tunnel junction (from Ref.[27j). Fluctuations 
of electric current, integrated over a short time interval At = 
5 ns and histogrammed, give the probability distribution of 
transmitted charge (insets). The mean, the variance and the 
skewness of this distribution are used to obtain the cumulants 
of FCS, Ci, C'2 and C3. The noise power S2, obtained from 
C2, agrees with the expectation for Poisson statistics, C2 = 
C\ (black line). The measured value of the 3rd cumulant is 
also close to the poissonian value, C3 ~ Ci (not shown, see 
Fig.3 in Ref. [2^]). Higher-order cumulants are challenging to 
obtain because the histogram of transmitted charge is nearly 
gaussian for nanosecond sampling times and currents of a few 
nA, as expected from the central limit theorem. 

Theoretical description of electron transport in a QPC 
involves scattering states constructed from the transmis- 
sion and reflection amplitudes A(t), B(t), which in gen- 
eral are time-dependent. In a Schrodinger representa- 
tion, the evolution of wavepackets undergoing scattering 
between the left and right reservoirs is described by 



U 



> 


L 




\x) 


R_ 





B(t x ) 
-A*{t x ) 



A(t x 
B*(t x 



\x(t))L 

\x{t)) R 



x < < x(t), 



and U\x)l,r. = \x)l,r otherwise. Here x(t) = x + vpt 
is the guiding center coordinate of a wavepacket, t x — 
—x/vf is the time of arrival at the scatterer, vp is the 
Fermi velocity, and \x}l,r describes incoming (x < 0) 
and outgoing (x > 0) wavepacket states in the leads. 

Crucially, the FCS generating function © can be ex- 
pressed through the evolved many-body density matrix 



of the initial state, projected on one of the leads after 
the evolution is completed, as in Eq.Q. This relation, 
derived in Refs. jl8ll34 |. is outlined below, and then used 
to find the entanglement entropy. 

There are certain general properties of the evolved den- 
sity matrix in a noninteracting fermion system which are 
best understood by considering evolution of a gaussian 
state 



■ exp 



X! // -.'"-".' 



Z = detfe""" + 1), (9) 



where H is a general hermitian operator in the single- 
particle Hilbcrt space, and Z is the normalization factor. 
The state of interest, describing the QPC at zero tem- 
perature, represents a particular case of Eq.©. 

One property which greatly simplifies the analysis, 
is that the state ([9]) remains gaussian under the evo- 
lution f4|. This follows from the observation that the 
Schrodinger evolution of ctj's is equivalent to the single- 
particle Heisenberg evolution of Hij. Gaussian form of 
the state is also preserved under projection [l8j |. 

This property can be used to reduce the many-body 
quantities of interest to certain one-particle quantities. 
Indeed, any gaussian state (J5J) can be described by a ma- 
trix in the single-particle Hilbert space defined as 



= Tr (pa\aA = 



1 



(10) 



In particular, for a fully filled Fermi sea in both reservoirs 
the matrix n is a projector on the subspace of all states 
with negative energy, el < 0, er < 0. As a projector, 
the matrix n satisfies the relation n 2 = n. 

In what follows, we will need to consider evolution of 
the projector n, followed by projection on the left reser- 
voir L, which is described by 



UnU\ M = P Lnu P L , 



(11) 



where Pp is a projection on the modes in L. The matrix 
njj is a projector describing evolved Fermi sea, whereas 
the matrix M, which is of main interest for us, is given 
by a product of three projectors. Thus generally M is 
not a projector. 

To illustrate the time evolution of single-particle quan- 
tities, such as njj, w e recall that in the FCS approach it 
is convenient to work in a time representation [43| . la- 
beling states by times of arrival at the scatterer t x . In 
this representation the initial Fermi projection is given 
by n(t,t') = 2iri(t-t'+io) I w ith / a 2 x 2 identity matrix 
in the L, R basis. The evolved state njj is given by 



nu(t,t') = U(t)n{t,t')U\t'), U{t) = 



Bit) A(t) 
-A*(t) B*(t) 



The evolution operator is diagonal with respect to the 
arrival time label t, which is precisely why this represen- 
tation is so convenient. 
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The FCS generating function ((6|) can be expressed in 
terms of the single-particle quantities, such as n, njj, Pl 
and M, in several different ways. The first representa- 
tion of this kind, found in Ref.[25|, involves a functional 
determinant 



X(A) = det(l -n + 



nU^e lXPL Ue 



i\p L 



)• (12) 



This determinant, which must be properly regularized for 
infinitely deep Fermi sea 31,; 3 2], can be explicitly eval- 
uated, yielding the FCS for various driving protocols of 
the QPC gHHi]. 

Another useful representation of the functional deter- 
minant giving x(A) was obtained recently in Ref.j34| (a 
similar relation was derived in a related problem of par- 
ticle number fluctuations (HHEI), 



X(A) = dct ( (1 - M + Me lX )e- lX{nPL)u 



(13) 



where (nP^u — UnPjJJ^ (see discussion in [3]). Here 
M is the Fermi sea in L, evolved by U and projected 
back to L by Pj, (see Eq. (fTT]) ). The unitary operator 
e -i\(nP L ) v con tributes a multiplicative factor of the form 
e lxX to the determinant, which may only affect the first 
cumulant C\ of the FCS generating function, Eq. ([8]) . 

The representation (fT3"|) is of interest because it reveals 
certain general features of x(A). It is convenient to in- 
troduce the spectral density of M, defined by [i(z) = 
Tr S(z — M). Since M is a product of three projec- 
tors, M = PhnuPh-, all its eigenvalues lie in the inter- 
val < z < 1 (see Refs. [33l39ll40| V Using the spectral 
density, we can rewrite (fT3"|) as 



log x(A) = ixX + I dzfx(z) log (l - z + ze iX ) . (14) 



This expression indicates that the FCS always assumes a 

generalized binomial form, x(A) EL (l ~ z + ze lX )^ z \ 
with the product taken over the entire spectrum of M. 



IV. THE RELATION BETWEEN 
ENTANGLEMENT AND FULL COUNTING 
STATISTICS 

The entanglement generated as a result of evolution of 
two Fermi seas coupled through a QPC is characterized 
by the von Neumann entropy of the density matrix ^ , 
S = — Tr pL{ti) logpi(ii). As we discussed above, this 
matrix is of a gaussian form, Eq. ([9]) . This property can be 
used to express entropy through single-particle quantities 
(see Refs. |18M4l| V giving 

S = -Tr (MlogM + (1 - M) log(l - M)) (15) 

where M is defined in (fTTj) , and the trace is taken in the 
space of single-particle modes in L. 



If the spectral density of M is known, the entropy fTS"]) 
can be written as 



S = - dzfi(z) (zlogz + (1 - z)log(l- z)) (16) 
Jo 

At the same time, as we discussed above, the FCS gen- 
erating function is also expressed through the spectral 
density of M . Furthermore, the spectral density is en- 
coded in the generating function. This can be seen most 
easily by rewriting Eq. (fT3"|) as 



X(z) = det ((z - M)e 



-i(nP L )u\(z) 



(1 



>)), (17) 



where the parameter A was changed to z = (1 — e ) . 
Thus the resolvent of M is given by the derivative 
d z log x(z — iO) up to a sum of two terms — + -^j arising 
from the last two factors in (fl~7|) . 

Using this observation the entanglement entropy can 
be expressed through the FCS cumulants C m [HI]. Writ- 
ing /j,(z) in Eq. lfTB)) as ^Im<9 z logx(z — *0) and plugging 
into it the relation ([8]), we integrate over z to obtain 



S 



m>0 



(2n) m \B r , 
0, 



m even 
m odd 



. (18) 



where B m are Bernoulli numbers, P>2 — ^, = — 
B 6 = j^... These numbers are defined by the generating 



function 



Z-~*r 



. B rl 



Interestingly, only even 



-1 ~ ^n=0 ' 

cumulants contribute to the entropy. The first few terms 
in the scries (IT? 



are: 



2 4 



2ir 6 



(19) 



2n! 



for large even n, which 
T8| stay bounded, 



Asymptotically, \B n \ ~ 
means that the coefficients in 
a m /m\ ~ 2 for large m. 

It is clear from this derivation that the relations (fTS]) 
and (|16p are completely general and valid for arbitrary 
driving. The formula (fl8|) can be used to determine the 
entanglement entropy from the values of FCS moments, 
whereas the formula (|16[) can be used when the spectral 
density /i(z) is known. 

As an example illustrating these relations, it is instruc- 
tive to consider a QPC biased with a DC voltage. The 
FCS for a DC-biased QPC with constant transmission 
< D < 1 is known to have binomial form (2f|, de- 
scribed by the generating function ([7]), 

X(A) = (l-D + De lX ) N , N = eVAt/h : (20) 

where N, given by the product of the bias voltage V and 
the measurement time At = t\ — to, is interpreted as the 
"number of attempts." Comparing this expression with 
(fTB")) . we infer that the spectral density of M in this case 
is a delta function, /i(z) = N6(z — D). Plugging this 



6 



result in the expression (fT6|) . we find that entanglement 
entropy is generated at a constant rate given by 

dS/dt = - (D log£> + (1 - D) log(l - D)) eV/h (21) 

Production of entanglement in a DC-biased QPC was 
considered in Ref . [42J , were the result (|2"Tj) was obtained. 
The rate of entanglement production is zero for D = 0, 1 
and maximal for D = 1/2. 

In this case, using the series (fT5)l turns out to be not 
too convenient because of a large number of high-order 
cumulants that contribute to the result. This can be 
seen most directly in the limit of small D <C 1, when the 
entanglement production rate scales as D log(l/D), while 
the cumulants C m ~ D. This means that there are about 
log(l/D) terms in the series (fl"8f giving contributions of 
the same order of magnitude. 

V. CONNECTING AND DISCONNECTING 
FERMI SEAS 

Here we shall discuss the protocol of driving the QPC 
in which it is opened at to < t < t\ and closed at t < to 
and ti < t, as illustrated in Fig[Tb. In the simplest case 
considered here, the QPC is unbiased, i.e. the right and 
left reservoirs remain at equal chemical potentials at all 
times. The fluctuations of charge transmitted through 
the QPC in this case are gaussian [i^ . 

x (A) = e-^, C 2 = i,log— , (22) 

TT Z T 

where At = t\ — to is the time window during which the 
QPC was open, and t is a short-time cutoff of order of 
the on/off switching time. 

The simplest way to estimate entanglement production 
is to use the formula (|18p. Since for a gaussian distribu- 
tion all cumulants are zero except C 2 , we find 

S = ^-C 2 = \\og—. (23) 
o or 

This result resembles the logarithmic dependence pre- 
dicted by conformal field theory, Eq.Q. In fact, since 
for free fermions the central charge is c = c = 1, the 
prefactor before the logarithm in Eq. §5§ is the same as in 
Eq.®. 

The similarity between the results (|23|) and ([3]) for en- 
tanglement entropy is of course not accidental. From a 
field-theoretic viewpoint, space and time play the same 
role in a conformal field theory. Therefore, analyzing en- 
tanglement using a window of size £ in space should be 
equivalent to doing it using a window of size At = l/vp in 
time, where Dp is Fermi velocity. This is precisely what 
the comparison of "time-like" Eq. (j2"3"|) and "space-like" 
Eq.© suggests. 

One can also understand the relation between the re- 
sults (I23[) and §3§ in a more intuitive way, without relying 
on a space-time duality. For that, we consider in more 



detail the process of mixing of two Fermi seas shown in 
FigHJa. Using different colors (blue and red) to mark par- 
ticles in different reservoirs, we observe that after t = t2, 
when the reservoirs are disconnected, there is a group of 
blue particles in the red Fermi sea. Simultaneously, there 
is a group of red particles in the blue Fermi sea. 

Assuming, without loss of generality, ballistic dynam- 
ics in each of the leads with constant velocity vf, we find 
that the blue and red groups of particles occupy spatial 
regions of size £ = VpAt. Since there is no correlation 
between the left and right reservoirs in the initial state 
of the system, we conclude that in the final state, shown 
in FigfJl there is no correlation between blue and red 
particles either, even if they reside in the same reservoir. 
This means that the von Neumann entropy evaluated for 
one of the reservoirs, in which both red and blue parti- 
cles are present simultaneously, will be the same as the 
entanglement entropy found for a single Fermi sea with 
a window of size £ = VpAt. 
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FIG. 4: Entanglement of two Fermi seas generated by con- 
necting and disconnecting them via a QPC during time in- 
terval At can be interpreted as entanglement in a stationary 
Fermi sea probed with a window of size £ = vpAt. 

This argument can in fact be made rigorous using the 
method of Ref.[43[. In this paper, concerned with the 
FCS of the process shown in Fig[2}D, the problem of fluc- 
tuations of charge transmitted through the QPC during 
a time interval At was mapped on the problem of particle 
number fluctuations in an interval of size £ = vpAt. Us- 
ing the bosonization representation of a free Fermi gas, 
the latter fluctuations can be shown to be gaussian, giv- 
ing Eq.(|32]), 

The result (|2"3"|) can be readily generalized to more com- 
plicated protocols of driving the QPC. In particular, it 
interesting to consider the QPC switching between the 

on and off states multiple times < t x < ... < t^ < 
t\ N \ In this case, generalizing the above argument, wc 
find gaussian charge statistics 

X (\)=e-^ c \ C 2 = -^G, (24) 

27T Z 

i,j=l L Q L l \ l l 

Using the relation (|18[) with the only nonvanishing con- 
tribution due to C2, we obtain the entropy 

7T 2 1 

S = —C 2 = -G. (26) 
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The case of multiple switching provides a time-like real- 
ization of the situation studied in Ref.[8J], where entan- 
glement of a conformal held theory was analyzed using 
a system of non-overlapping windows < x^ 



< 



< 



x$ < x[" ' . As in the case of single switching, Eq.Q, 





the answers for the entanglement entropy coincide after 

M 
'o.i 



v F t 01 . 



distances are converted to times via x 

Multiple switching of the QPC between the on and off 
states, repeated periodically in time, can be used to con- 
vert the current fluctuations, described by Eqs. (f2"2"l) . (|2"4")) . 
into the DC shot noise which can be measured by conven- 
tional techniques. In Ref.[l8[, the effective temperature 
of such noise was estimated to be about 25 mK for the 
driving frequency v = 500 MHz, putting it in the experi- 
mentally feasible range. 

In view of the possibility of such an experiment, it is 
interesting to consider how the above results are changed 
if the QPC transmission in the open state is less than one. 
This problem can be readily addressed both for single and 
multiple switching protocols, since the FCS in this case 
has been known [431 ] . 

X (A) = e-^ G , sin ^A* = y/D sin ~ A, (27) 



with G given by ([25]) as above, and D < 1 the QPC 
transmission coefficient. 

Notably, in this case x(A) is non-gaussian. Thus the 
simplest way to find the entropy is to use its relation 
with the spectral density of M, Eq. ([T6[) . The spec- 
tral density fj,(z) can be evaluated using the resolvent 
Tr (l/(z - M - i0)), which is found from x(A), Eq.([27]), 
as discussed above. This gives a peculiar function [18| 
which vanishes inside the interval 

f 



Z\ < Z < 2 2 , ^1,2 = - \1 T VT-Dj , 

while outside this interval, at < z < z\ and z<i < z < 1, 
it is given by 



n(z) 



G 



II - 2z\VD 



^Z{l-Z)y/{Z-Zl){z-Zi) 



(28) 



The entropy, found from (|16p . exhibits a logarithmic de- 

pendence on the times t \ identical to (|25j) ; the only 
change is a .D-dependent prefactor. 

The behavior of the factor F — S(D)/S(1) describ- 
ing entropy reduction due to imperfect transmission in 
the QPC was analyzed in Ref . [l8j . It was found that, 
unless D is very small, the reduction in entropy can be 
attributed mostly to the change in the second cumulant, 
C*2 = 2^^' w ith the contribution from higher cumu- 
lants being relatively small. Thus even for imperfect 
QPC transmission, the DC shot noise generated by QPC 
switching gives a reasonably good estimate of the entan- 
glement entropy production. 



VI. CONCLUSIONS 

There are several ways in which the relation between 
many-body entanglement and the Full Counting Statis- 
tics (FCS) description of quantum transport comes as 
noteworthy and surprising. First, it provides an inter- 
esting new application for the FCS approach. Using the 
framework outlined above, many results from the FCS 
literature can be reinterpreted and used to study entan- 
glement production in a variety of regimes of experimen- 
tal interest. 

Second, the FCS approach offers new insight into the 
nature of many-body entanglement in driven quantum 
systems. In particular, two different kinds of entan- 
glement, called restricted and unrestricted, are distin- 
guished in quantum information theory [42J. This dis- 
tinction refers to our ability to detect entanglement by 
means of local measurements, as required in various tele- 
portation and cryptography protocols. Entanglement 
production in a system such as the QPC, opened and 
closed multiple times as discussed above, results from 
particle exchange between two Fermi seas. However, 
since particle number is a conserved quantity, some of 
the produced entanglement may be inaccessible to local 
measurements, and thus unuseful from a quantum infor- 
mation standpoint. 

How large is the unrestricted ( "useful" ) entanglement? 
This question can be addressed by generalizing the FCS- 
based approach, as discussed in detail in Ref. [44j . It 
turns out that for realistic driving protocols, such as 
those analyzed above, nearly all entanglement is unre- 
stricted. Specifically, when the number of particles trans- 
ferred through the QPC becomes large, the restricted en- 
tanglement entropy scales as log S, where S ^> 1 is the 
total entanglement entropy. Thus only a small fraction 
of entanglement is degraded to the restricted form due 
to particle conservation. 

Finally, the relation between entanglement and FCS 
opens a way to perform direct measurement of entangle- 
ment entropy by detecting fluctuations of electric cur- 
rent in a driven system. In particular, by using the QPC 
switching periodically between the on and off states, and 
utilizing space-time duality of one-dimensional systems, 
the relation between entanglement production and noise 
can be used to test the seminal S = | log L prediction of 
conformal field theory More generally, this relation 
offers a method for theoretical and experimental investi- 
gation into the nature of many-body entanglement, and 
in particular, of its build up in non-equilibrium quantum 
systems. 

We thank C. W. J. Beenakker, B. Reulet, and M. 
Reznikov for useful discussions. This research was sup- 
ported in part by US-Israel Binational Science Founda- 
tion and W. M. Keck Foundation Center for Extreme 
Quantum Information Theory (L.L.). 
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